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Abstract 


We report the first measurement of the absolute branching fraction for At — Ани. This measurement is based on 
a sample of e* e^ annihilation data at a center-of-mass energy of \/s = 4.6 GeV collected with the BESIII detector at 
the BEPCII storage rings. The sample corresponds to an integrated luminosity of 567 pb~!. The branching fraction 
is determined to be B(A* > Au*v,) = (3.49 + 0.46(stat) + 0.27(syst))%. In addition, we calculate the ratio 
B(A} > Aptv,)/B(AE > Aetve) to be 0.96 + 0.16(stat) + 0.04(syst). 


Keywords: At, semi-leptonic decay, absolute branching fraction, BESIN 


1. Introduction 


Semileptonic (SL) decays of the lightest charmed 
baryon, AZ, provide a stringent test for non-perturbative 
aspects of the strong interaction theory. The At — 
Лету, (L denotes lepton) decay is dominated by the 
Cabibbo-favored transition с — s£* v;, which occurs 
independently of the spin-zero and isospin-zero spec- 
tator ud diquark, to good approximation. This leads 
to a simpler theoretical description and greater predic- 
tive power in the non-perturbative models than in the 
case for charmed mesons [1]. Predictions of the branch- 
ing fraction (BF) B(A} — A€tve) in different the- 
oretical models vary over a wide range from 1.4% to 
9.296 [2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13], depending on 
the choice of A^ wave function model and the treatment 
of decay dynamics. In 2015, BESIII measured the ab- 
solute BF for A^ — Лети, to be B(A} > Aetve) = 
(3.63 + 0.38 + 0.20)% [14], which disfavors the predic- 
tions in Refs. [2, 3, 5, 6, 7] at the 9596 confidence level. 
It is desirable to test these theoretical predictions by 
measuring B(A} — Лито). In addition, lepton uni- 
versality can be tested by comparing the BFs between 
the electronic and muonic modes. 

In this paper, we report the first absolute measure- 
ment of B(A} — Ag*v,,) by analyzing a data sample 
with an integrated luminosity of 567 pb~! [15] accu- 
mulated at a center-of-mass (c.m.) energy of ys = 
4.6 GeV with the BESIII detector at the BEPCII col- 
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lider, which is the largest етет collision sample near 
the A+ A> mass threshold. At this energy, the At is 
produced in company with one A> baryon only, and 
no other hadrons are kinematically allowed. Hence, 
B(A} — Aptv,) can be accessed by measuring the 
relative probability of finding the SL decay when the 
A; is detected in a number of prolific decay channels. 
This will provide a straightforward and direct BF mea- 
surement without requiring knowledge of the total num- 
ber of A+ A> pairs produced. In the following, charge 
conjugated modes are always implied, unless explicitly 
mentioned. 


2. BESIII Detector and Monte Carlo Simulation 


The BESIII [16] detector is a cylindrical detector 
with a solid-angle coverage of 93% of 4r that oper- 
ates at the BEPCII collider. It consists of a Helium- 
gas based main drift chamber (MDC), a plastic scintilla- 
tor time-of-flight (TOF) system, a CsI (ТІ) electromag- 
netic calorimeter (EMC), a superconducting solenoid 
providing a 1.0 T magnetic field and a muon counter. 
The charged particle momentum resolution is 0.5% at a 
transverse momentum of 1 GeV/c. The photon energy 
resolution in the EMC is 2.5% in the barrel and 5.0% 
in the end-caps at 1 GeV. More details about the design 
and performance of the detector are given in Ref. [16]. 

A GEANTA-based [17] Monte Carlo (MC) simu- 
lation package, which includes the geometric descrip- 
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tion of the detector and the detector response, is used 
to determine the detection efficiency and to estimate 
the potential backgrounds. Signal MC samples of a 
АТ baryon decaying only to Aptv, together with a 
Az decaying to specified modes are generated with the 
KKMC [18] and EVTGEN [19], taking into account the 
initial state radiation (ISR) [20] and the final state ra- 
diation (FSR) [21] effects. For the simulation of the 
process At — Ag*v,, we use the form factor ob- 
tained using Heavy Quark Effective Theory and QCD 
sum rules in Ref. [10]. To study backgrounds, inclu- 
sive MC samples are simulated which consist of At A7 


events, DD + X production, ISR return to the 
charmonium(-like) ~ states at lower masses, and QED 
processes. The decay modes with known BFs of the 
Ac, Y and D, particles taken from Particle Data Group 
(PDG) [22] are simulated with EVTGEN, while the re- 
maining unknown decays are generated with LUND- 


CHARM [23]. 


3. Analysis 


Following the similar technique of the single tag 
(ST) and double tag (DT) in Ref. [14], we reconstruct 
the A> baryons in eleven hadronic decay modes as 
listed in the first column of Table 1. The intermedi- 
ate particles КО, A, X09, У- and л? are reconstructed 
through their decays КО — път, A > prt, 309 > 
ҮЛ with A > pat, X- — pr? and ло — yy, re- 
spectively. The detailed selection criteria for charged 
and neutral tracks, 7°, КО and A candidates used in the 
reconstruction of tags are described in Ref. [14]. 

In this analysis, the ST Л; signals are identified 
using the beam energy constrained mass, Mpc = 


"ES а/с — |0 - |2/е2, where Epeam is the beam 


energy and PA- is the momentum of the A7 candi- 
date. To improve the signal purity, the energy differ- 
ence AE = Ebpbeam — Аг for each candidate is re- 
quired to be within +30, g around the AF peak, where 
CAE is ће AF resolution and E Аг is the reconstructed 


Л7 energy. Table 1 shows the mode dependent AF re- 
quirements and the ST yields in the Mgc signal region 
(2.280, 2.296) GeV/c”, which are obtained by a fit to the 
Mec distributions. The detailed process to extract the 
ST signal yields is described in Ref. [14]. The total ST 
yield summed over all 11 modes is N$% = 14415+159, 
where the uncertainty is statistical only. 

Candidate events for A^ — Ag, are selected 
from the remaining tracks recoiling against the ST A7 
candidates. The A candidate is formed from a рт com- 
bination that is constrained by a common vertex fit to 


have a positive decay length L. If multiple A candidates 
are formed, the one with the largest L/or, is retained, 
where or, is the resolution of the measured L. Parti- 
cle identification (PID) is performed using probabilities 
derived from the specific energy loss dE/dx measured 
by the MDC, the time of flight measured by the TOF, 
and energy measured by the EMC; a и candidate is re- 
quired to satisfy C, > 0.001, Li, > £z and Li, > Ly, 
where Lis Li, and £^. are the probabilities for a muon, 
electron, and kaon, respectively. 

Studies on the inclusive MC samples show that the 
backgrounds are dominated by At — Ant, X?z* and 
Ат+т°. Backgrounds from At — Amt and At > 
У0т+ are rejected by requiring the Ли? invariant mass 
Мл „+ is less than 2.12 GeV/c?. The background from 
At — Artn” is suppressed by requiring the largest 
energy of any unused photons Emax be less than 0.25 
GeV and the deposited energy for the muon candidate 
in the EMC be less than 0.30 GeV. 

Since the neutrino is not detected, we employ the 
kinematic variable Umiss = Emi — |Dmiss c| to iden- 
tify the neutrino signal, where Ei, and [мв аге 
the missing energy and momentum carried by the neu- 
trino, respectively. They are calculated as Emiss = 
Exeam — Ёл = E, and Pmiss = Pat = DA = Pats 
where p,+ is the momentum of the Af baryon, EA(DA) 
and E,,+(p,,+) are the energies (momenta) of the A and 
L^, respectively. Here, the momentum р A; is given 


by Dat = —Ptag ае 


momentum direction of the ST A> and ту is the nom- 


2 EE 
=mi where tag is the 


inal AI mass [22]. For the signal events, the Umiss dis- 
tribution is expected to peak at zero. 

The distribution of the рт invariant mass M,,,- 
versus Umiss for the A^ — Aptv, candidates in data 
is shown in Fig. 1 (a), where a cluster around the signal 
region is evident. After requiring M,,,- to be within 
the A signal region, the projection of Umiss is shown in 
Fig. 1(b). Two bumps, which correspond to the signal 
peak (left side) and background At — Az*7? (right 
side), are visible. According to MC simulations, the 
survival rate of the background process At — Az a? 
is estimated to be 7лл+ хо = (3.67 + 0.05)%, where the 
BFs for A + рл and 1? — yy are included. Thus, the 
number of ће A} — Az? background events can be 
estimated by 


МЕК о = Ny B(AF > Artn?) Mantan: (1) 


Inserting the values of NP, ПАлт+хо and B(A} — 
Artn?) = (7.01 + 0.42)% [24] in Eq. (1), we obtain 
NPE a = 37.1 + 2.3, where the uncertainties from 


TABLE 1: AF requirements and ST yields N- АТ in data, in which the uncertainties are statistical only. 


Mode AE (GeV) Ni 
PKs [—0.025, 0.028] 1066 + 33 
рК т^ [—0.019, 0.023] 5692 + 88 
РКО? [—0.035, 0.049] 593 + 41 
pK^nz-4? [—0.044, 0.052] 1547 + 61 
pKintn- [-0.029,0.032] 516+ 34 
Ат [—0.033, 0.035] 598 + 25 
Ат п0 [—0.037, 0.052] 1864 + 56 
Arnt [—0.028, 0.030] 674+ 36 
Lone [—0.029, 0.032] 532 + 30 
Er? [—0.038, 0.062] 329 + 28 
X-mtm- [—0.049, 0.054] 1009 + 57 
Nx, Nan+no and B(At — An*7®) are included. with 
We apply a fit to the Umiss distribution to obtain the x Р NN ‚ 
signal yields. The A] — Лито, signal shape is de- B(A? > Ашти) = NES pu, BAS gr) 
Ac Бег 


scribed with a function f, which consists of a Gaussian 
function to model the core of the Umiss distribution and 
two power law tails to account for the effects of ISR and 
FSR in the form [25] of 


pac ai Tt)", to 
f (Umiss) = ef 2/2, —ag < t < ay. 
TOE 2, t<-ap 


(2) 
Here, t= (Umiss = Umean)/OUmiss> Umean and OU riss 
are the mean value and resolution of the Gaussian func- 
tion, respectively, pı = (nı/a1)”1e721/2 and рә = 
(na/a3)n2e-93/2, The parameters a1, 2, n1 and n 
are fixed to the values obtained by studying the signal 
MC simulations. For backgrounds, a double Gaussian 
function with parameters fixed according to MC simu- 
lations is used to describe the A^ — Az * x? peaking 
background and a MC-derived shape is used to describe 
other combinatorial backgrounds. In the fit, we fix the 
number of ће A+ — Az? background events to be 
estimated №» d „о as described above. From the fit, we 
obtain the ае of events of АТ — Aptv, to be 


NV, = = 78.7 + 10.5, where the киш is sta- 


tistical only. A fit with unconstrained NP ДЕ: Ж gives 
77.1 + 11.4 events of signal, which is in good agree- 
ment with the estimation when NES a is fixed. Based 
on the data in A sidebands in Fig. 1(a), the background 
events from the non-A SL decays is found to be negligi- 
ble. 


The absolute BF for АТ — Aptv, is determined 


(3) 
where €,,,+,,, is the detection efficiency for the At > 
Aptv, decay, which does not include the BF for A > 
pm. For each ST mode 7, the efficiency E jin, is ob- 
by the 
ST efficiency Eur After weighting є* Au Dy with the ST 
yields in data for each ST mode 7, we determine the 
overall average efficiency &4,,«,, = (24.5 + 0.2)%. 
By inserting the values of aus з М, EAutv, and 
B(A — рт) [22] in Eq. (3), we obtain B(A} > 
Лито) = (3.49 + 0.46 + 0.27)%, where the first un- 
certainty is statistical, and the second uncertainty is sys- 
tematic as described below. 

With the DT technique, the uncertainties on the BF 
measurement are insensitive to those originating from 
the ST side. The systematic uncertainties for measur- 
ing B(A} — Aptv,) mainly arise from the uncertain- 
ties related to the tracking and PID of the muon can- 
didate, A reconstruction, Umiss fit, peaking background 
subtraction, Fymax and M,,,+ requirements, and sig- 
nal MC modeling. Throughout this paragraph, the sys- 
tematic uncertainties quoted are relative uncertainties. 
The uncertainties of the џ+ tracking and PID are deter- 
mined to be 1.0% and 2.0%, respectively, by studying 
a control sample of e*e^ — (y)u*p~ events. The 
uncertainty of the A reconstruction is determined to be 
2.5% by studying a control sample of xez — AArt r7 
decays. The uncertainty of Umiss fit is estimated to 
be 1.5% obtained by varying the fitting range and ex- 
amining the fluctuation of the non-peaking background 
shape. The uncertainty due to peaking background 


tained by dividing the DT efficiency figi Хи 
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FIG. 1: (a) Distribution of Mpr- versus Umiss for the At — Aptv, candidates in data. The area between the dashed lines denotes the A 
signal region and the hatched areas indicate the A sideband regions. (b) Fit to the Umiss distribution within the A signal region. Data are shown as 
the dots with error bars. The long-dashed curve (green) shows the At — Antr’ background while the dot-dashed curve (blue) shows other At 


decay backgrounds. The thick line (red) shows the total fit. 


At — Amm? subtraction is estimated to be 2.5% ob- 
tained by varying Ne „o equivalent variations of +10 
of the quoted BFs, and smearing the MC-derived shape 
of At — Artnr’ backgrounds with a Gaussian func- 
tion to accommodate the resolution difference between 
the data and MC simulation. The uncertainty in the 
Eymax requirement is estimated to be 2.6% by using a 
control sample of e* e^ — ppt events. The uncer- 
tainty in the МА „+ requirement is estimated to be 2.0% 
by comparing the obtained B(A} — Aptv,,) under ће 
alternative requirements of M4,« < 2.07 GeV/c? or 
Млн+ < 2.17 GeV/c? with the nominal value. The 
uncertainty in the signal MC model is estimated to be 
5.2% by varying the parameterization of the form factor 
function according to Refs. [10, 26] and by taking into 
account the q? dependence observed in data. In addi- 
tion, there are systematic uncertainties from the quoted 
B(A — pr7) (0.8%), the Ме (1.0%) evaluated by us- 
ing alternative signal shapes in the fits to the Мвс spec- 
tra [14], and MC statistics (0.8%). All these systematic 
uncertainties are summarized in Table 2, and the total 
systematic uncertainty is evaluated to be 7.7% by sum- 
ming up all the individual contributions in quadrature. 


The ratio of branching fractions B(A} — 
Apty,)/B(At ә Aetve) is calculated with the 
measured B(A} —› Apty,) in this work and 
B(At — Aety,) = (3.6340.38(stat) +£0.20(syst))% 
from ВЕЅШ [14]. We determine B(A} — 
Aptyv,)/B(At — Ле+и,) to be 0.96 + 0.16 + 0.04, 
where the first uncertainty is statistical and the second is 
systematic, in which the common systematic uncertain- 
ties from the tracking efficiency, the A reconstruction, 
the quoted BF for A — ртт, the number of A; tags 
NS and the MC model cancel. 


4. Summary 


In summary, based on the e* e^ collision data corre- 
sponding to an integrated luminosity of 567 рь! taken 
at \/5 = 4.6 GeV with the BESIII detector, we re- 
port the first direct measurement of the absolute BF for 
AT + Аии, to be (3.49 + 0.46 + 0.27)%, where the 
first uncertainty is statistical and the second is system- 
atic. The result is consistent with the value in PDG [22] 
within 2c of uncertainty, but with improved precision. 
This study helps to extend our understanding on the de- 
cay mechanism of the A7 SL decay. Based on this 
result and the previous BESIII work [14], we deter- 
mine the ratio B(A} > Aptvy,)/B(At > Ает) = 
0.96 + 0.16 + 0.04, which is compatible with unity. As 
the theoretical predictions on B(A} — A£*v;) vary in 
a large range of 1.4% to 9.2% [2, 3, 4, 5, 6, 7, 8, 9, 
10, 11, 12, 13], the measured B(A? — Aptv,,) in this 
work апа B(A} — Ae*y,) in Ref. [14] provide strin- 
gent tests on these non-perturbative models. 
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